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Abstract. A general method for treating non-Gaussian wave functionals in quantum field theory is presented and applied to 
the Hamiltonian approach to Yang-Mills theory in Coulomb gauge in order to include a three-gluon kernel in the exponential 
of the vacuum wave functional. The three-gluon vertex is calculated using the propagators found in the variational approach 
with a Gaussian trial wave functional as input. 
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INTRODUCTION 

Over the last few years, there have been substantial ef- 
forts devoted to a variational solution of the Yang-Mills 
Schrodinger equation in Coulomb gauge [1, 2, 3]. In this 
approach, using Gaussian type wave functionals, mini- 
mization of the energy density results in the so-called 
gap equation for the inverse equal-time gluon propagator. 
This equation has been solved analytically in the ultravi- 
olet [2] and in the infrared [4], and numerically in the 
full momentum regime [2, 3]. One finds an inverse gluon 
propagator which in the UV behaves like the photon en- 
ergy but diverges in the IR, signalling confinement. The 
obtained propagator also compares favourably with the 
available lattice data. There are, however, deviations in 
the mid-momentum regime (and minor ones in the UV) 
which can be attributed to the missing gluon loop, which 
escapes the Gaussian wave functionals. These deviations 
are presumably irrelevant for the confinement properties, 
which are dominated by the ghost loop (which is fully in- 
cluded under the Gaussian ansatz), but are believed to be 
important for a correct description of spontaneous break- 
ing of chiral symmetry [5]. 

In this talk, we present a generalization of the varia- 
tional approach to the Hamiltonian formulation of Yang- 
Mills theory [2] to non-Gaussian wave functionals. The 
expectation value of the Hamilton operator can be ex- 
pressed in terms of the variational kernels occurring in 
the ansatz through Dyson-Schwinger equations (DSEs). 
The three-gluon vertex and the effects of the gluon loop 
on the gluon propagator are investigated. 

NON-GAUSSIAN WAVE FUNCTIONALS 

In the Hamiltonian approach to Yang-Mills theory in 
Coulomb gauge, the vacuum expectation value (VEV) 



of an operator depending on the transverse gauge field 
is given by 

(K[A})= J mjAlvlA^KlA], (1) 

where J a = Det(G,4) is the Faddeev-Popov determinant, 
Ga = (—dD)^ 1 is the inverse Faddeev-Popov operator, 
and y/[A] is the vacuum wave functional. In Eq. (1), the 
functional integration runs over transverse field configu- 
rations djA" = and is restricted to the first Gribov re- 
gion. Writing the vacuum functional as 

|V/[A]| 2 =:exp{-S[A]}, (2) 

one can derive DSEs from the functional identity 

0=J^A{7 A e- s %[A]}. (3) 

In this talk, we consider a functional of the form 

S{A] = coA 2 + ^y 3 A\ (4) 

where a> and 73 are variational functions. A vacuum 
functional containing also a quartic term is discussed 
in Ref. [6]. With the explicit form Eq. (4) for the vac- 
uum functional, the DSEs are the usual DSEs of Landau 
gauge Yang-Mills theory, however, in D = 3 dimensions 
and with the bare vertices of the usual Yang-Mills ac- 
tion replaced by the variational kernels. For the gluon 
(AA) =: 1 /(2£2) and ghost propagator (Ga) these equa- 
tions are shown in Fig. 1. It should be stressed that these 
Hamiltonian DSEs are not equations of motion in the 
usual sense, but rather relations between the Green func- 
tions and the so far undetermined variational kernels. 
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FIGURE 1. Dyson-Schwinger equation for the gluon (top) 
and ghost propagator (bottom). Here and in the following, 
small filled dots represent propagators, small empty dots vertex 
functions, and empty boxes the variational kernels. 



VARIATIONAL APPROACH 

The Yang-Mills Hamilton operator in Coulomb gauge 
reads [7] 



r r i _j s s i 2 

H ™-J ~2 Ja SA Ja SA + 2 B 



(5) 



where B is the non-abelian magnetic field, A is the gauge 
field in the adjoint representation of the colour group, and 
Fa = Ga (— <5 2 )Ga is the Coulomb interaction kernel. The 
vacuum energy is evaluated as VEV of the Hamilton op- 
erator Eq. (5) with the vacuum state defined by Eqs. (2) 
and (4). By using the DSEs stemming from the identity 
Eq. (3), the energy density can be written as a functional 
of the variational kernels, 



(H YM ) =E[a>,y i \. 



(6) 



By using a skeleton expansion, the vacuum energy can be 
written at the desired order of loops. Confining ourselves 
to two loops, the variation of the vacuum energy Eq. (6) 
with respect to the kernel 73 fixes the latter to 



2igf abc 

Sij{p-q)k + 8 jk (q- 



k)i + 8 ki (k-p)j 



£!(£)+ £1(5) + £!(*) 



(7) 



Equation (7) is reminiscent of the lowest-order pertur- 
bative result [8], with the perturbative gluon energy \p\ 
replaced by the non-perturbative one Cl(p). 

Combining the gluon DSE with the variational equa- 
tion for the two-gluon kernel CO, one arrives at the gap 
equation for the gluon propagator 



n(py=p 2 +x(pr+ic(p)-iG(p). 



(8) 



The explicit expressions for the loop terms can be found 
in Ref. [6]. With a Gaussian wave functional, only 
the ghost loop %{p) and the contribution Ic(p) of the 
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FIGURE 2. Gluon propagator obtained with a Gaussian 
(dashed line) and a non-Gaussian functional (straight line), 
compared to the lattice data from Ref. [9]. 
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FIGURE 3. Truncated DSE for the three-gluon vertex, under 
the assumption of ghost dominance. 



Coulomb kernel appear in the gap equation. The pres- 
ence of the gluon loop Ia(p) in the gap equation (8) mod- 
ifies the UV behaviour and allows us to extract, from the 
non-renormalization of the ghost-gluon vertex, the cor- 
rect first coefficient of the j3 function. The solution of 
the full set of coupled integral equations for the ghost and 
gluon propagators and for the Coulomb form factor (see 
Ref. [6]) is in progress. Here, in order to estimate the size 
of the gluon-loop contribution to the gluon propagator, 
we use the gluon and ghost propagators obtained with a 
Gaussian wave functional [3] to calculate the gluon loop. 
The result is shown in Fig. 2, together with lattice data 
from Ref. [9]. The agreement between the continuum 
and the lattice results is improved in the mid-momentum 
regime by the inclusion of the gluon loop, i.e. the three- 
gluon vertex, as observed also in Landau gauge [10]. The 
mismatch in the UV is a consequence of the approxima- 
tions involved, and should disappear when the full sys- 
tem of coupled equations is solved. 



THREE-GLUON VERTEX 

The truncated DSE for the three-gluon vertex under 
the assumption of ghost dominance is represented dia- 
grammatically in Fig. 3. Possible tensor decompositions 
of the three-gluon vertex are given in Ref. [11]. Here, 
for sake of illustration, we confine ourselves to the form 
factor corresponding to the tensor structure of the bare 
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FIGURE 4. Form factor f^A of the three-gluon vertex for 
orthogonal momenta and comparison to lattice data for the d = 
3 Landau-gauge vertex [12]. The momentum scale is arbitrary 
and has been adjusted to make the sign change occur at the 
same point. The lattice data are shown by courtesy of A. Maas. 

three-gluon vertex 

ryrff 

/3A: ~ r (o) r (o)' (9) 

where rl ^ is the perturbative vertex, given by Eq. (7) 
with £l(p) replaced by \p\. Furthermore, we consider a 
particular kinematic configuration, where two momenta 
have the same magnitude 

P 2 l=P2=P 2 , Pl-P2 = C P 2 . (10) 

To evaluate the form factor /t,a(p 2 ,c), we use the ghost 
and gluon propagators obtained with a Gaussian wave 
functional [3] as input. The IR analysis of the equation 
for /3a(p 2 ,c) [Eq. (9)] performed in Ref. [4] shows that 
this form factor should behave as a power law in the IR, 
with an exponent three times the one of the ghost dress- 
ing function; this is confirmed by our numerical solution 
[6]. The result for the scalar form factor f^A for orthogo- 
nal momenta, f(p 2 ,0), is shown in Fig. 4, together with 
lattice results for d = 3 Landau gauge Yang-Mills the- 
ory. Our result and the lattice data compare favourably 
in the low-momentum regime. In particular, in both stud- 
ies, the sign change of the form factor occurs roughly at 
the same momentum where the gluon propagator has its 
maximum. (The scale in Fig. 4 is arbitrary.) 

CONCLUSIONS 

We have presented a method to treat non-Gaussian wave 
functionals in the Hamiltonian formulation of quantum 
field theory. By means of Dyson-Schwinger techniques, 



the expectation value of the Hamiltonian is expressed in 
terms of kernels occurring in the exponent of the vac- 
uum wave functional. These kernels are then determined 
by minimizing the vacuum energy density. We have esti- 
mated the three-gluon vertex by using the propagators 
found with a Gaussian wave functional as input. The 
result compares fairly well to the available lattice data 
obtained in d = 3 Landau gauge. The gap equation for 
the gluon propagator contains the gluon loop, which was 
missed in previous variational approaches with a Gaus- 
sian wave functional. The gluon loop gives a substantial 
contribution in the mid-momentum regime while leaving 
the IR sector unchanged, and it also provides the cor- 
rect asymptotic UV behaviour of the gluon propagator 
in accord with perturbation theory [8]. The presently de- 
veloped approach allows a systematic treatment of cor- 
relators in the Hamiltonian formulation of a field theory, 
and opens up a wide range of applications. In particular, 
it allows us to extend the variational approach from pure 
Yang-Mills theory to full QCD. 
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